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1 Introduction 

"Charming penguins" started back in 1997, coming out of a study aimed to evaluate 
hadronic effects in non-leptonic two-body B decays. During the years, several episodes 
added to the saga: 

1. The penguin menace [[[], |2|. A lattice- inspired Wick-contraction parametrization 
of hadronic amplitudes was introduced and the observation was put forward that 
non-factorizable penguin contractions of current-current operators containing 
two c quarks (the charming penguins) could give large contributions in some B 
decay channels, notably B — > Kn (a similar idea was already present in ref. 0). 

2. The neat hack of the clones [|J. The original Wick-contraction parametrization 
was modified by Buras and Silvestrini. The hadronic matrix elements were ex- 
pressed in terms of new renormalizat ion-group invariant parameters given by 
suitable combinations of the old ones. Many B decay channels were classified 
according to the new parametrization. Charming penguins became a more com- 
plex object, containing further contractions (annihilations, penguin contractions 
of penguin operators) in addition to the original one. 

3. A new hope []5], |6|]. The one-loop proof that factorization of hadronic matrix 
elements holds in the limit — > 00 puts phenomenological approaches based 
on factorization on a firmer theoretical ground (other theoretical approaches 
to factorization in the infinite mass limit were already developed, although 
not at the same level of accuracy 0). In this limit, non-factorizable corrections 
were shown to be computable using perturbation theory. Perturbative penguins 
turned out to give in general small contributions. Charming penguins seemed 
at loss. 



1 Talk given by M. Ciuchini. 



4. Charming penguins strike back |§. Using B — > Ktt data, it was shown that 
the parameter accounting for charming penguins has the expected size of a 
h-QCD/™>b correction. Therefore, a sizable non-perturbative effect of charming 
penguins is not in disagreement with the results on factorization obtained in 
the infinite mass limit. In addition, it is preferred by the data. 

5. The return of factorization H[K|. While everybody agrees that power-suppressed 
terms are in general non-perturbative and non-factorizable, it was argued that 
still the bulk of the Agc^/m?, corrections can be either factorized or, failing 
that, accounted for by few parameters (this framework is called improved QCD 
factorization). In addition, these parameters, once properly defined, are claimed 
to have negligible effects on B — > Ktt branching ratios. Under these assump- 
tions, which were shown to be compatible with the present data, these branching 
ratios can be used to extract the CKM angle 7. 

Is the saga arrived to its end? Theoretically, it is not clear whether a non-perturbative 
contribution such as charming penguins is large or small. A recent calculation using 



renormalons found no sign of it |Tl| , while, on the contrary, it is present and effective 
in other approaches 1 12| . The B — > Ktt data certainly call for power-suppressed 
terms and charming penguins are able to provide what is needed. Other approaches, 
such as the popular improved QCD factorization, are also compatible with the data, 
but none is able to make predictions based only on the theory, due to the presence of 
phenomenological parameters. Indeed, the presence of these parameters makes us very 
skeptical about the possibility of extracting the CKM angle 7 from the measurement 
of the B — > Ktt branching ratios. 



2 Charming penguins at work 

In this section we collect the main formulae for the amplitudes of B — > Ktt, tttt, 
introducing the parametrization used in the analysis. We refer the reader to the 
literature for any detail on the origin and the properties of these parameters |], |[ f|, §]. 
From ref. H, one reads 

A{B d ^K+TT-) = ^{KP 1 -K{E 1 -P^)) 

A(B + ^K + tt°) = ( ^(\ s t P 1 -\ s u {E l + E 2 -P^ IM + A l )) + AA 

A(B + ^K°tt + ) = ~^=( ~ K^i + K.(Ai - -Pi G/M )) + AA 

A(Bd — > K°TT ) = ^^-y tPl -X s u (E 2 + P^ M ))+AA (1) 
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where X^, = V q ' q V q * b . Neglecting the A iy these parameters can be rewritten as 

E 1 = a1A nK , E 2 = a c 2 A Kn , A x = A 2 = , 

Pi = a\A^ K + Pi , P X IM = [a\ — aX)A^K + P 1 G/M . 



A( B, k + x- ) -| (\ d t (Pi + P 3 ) - XHEk + A 2 -P[ 
A(B d — >■ 7r + 7T°) = ^(-A^(Pi + P 2 ))+AA 
A(P d ^7r°7r°) 



A f s (Pi + P 3 ) - A:(P 2 + Pf ZM + P 3 GIM - A 2 )) + AA , 



(2) 



The terms proportional to ai gives the parameters computed in the limit m& — > oo 
using QCD factorization. Their definition, together with those of A^ K , A Kw , etc., 
can be found for instance in ref. |J, although power-suppressed terms included there, 
proportional to the chiral factors should be discarded in eqs. (Q). In our case, in 
fact, terms of O^AqcD/^b) are accounted for by two phenomenological parameters: 
the charming-penguin parameter Pi and the GIM-penguin parameter p^ IM . In 
B — > Kix there are no other contributions, once flavour SU(2) symmetry is used 
and few other doubly Cabibbo-suppressed terms, including corrections to emission 
parametes Pi and E 2 , some annihilations (Ai) and the Zweig-suppressed contactions 
(AA), are neglected M. On the contrary, further power-suppressed terms (A 2 , P 3 , 
P^ IM ) enter the B — > tttt amplitudes, all with the same power of the Cabibbo angle. 
Therefore, these modes are subject to a larger uncertainty than the B — > Kt\ ones. 

Using the inputs collected in Table [l], we fit the value of the complex parameter 
Pi = (0.13 ± 0.02) e^ 114 * 35 )" in units of / 7r P 7r (M^). Notice that the sign of the phase 
is practically not constrained by the data. This result is almost independent of the 
inputs used for the CKM parameters p and 77, namely whether these parameters are 
taken from the usual unitarity triangle analysis (UTA) [|14], |TS[ or only the constraint 
from iKift/Vc&l is used. 



KblxlO 3 



KJxio 3 



B K 



f Bd ^B~ d (MeV) 



40.9±1.0 

F K {Ml\ 



3.70±0.42 

B(K+tt-)x10 6 



0.86±0.06±0.14 230±30±15 
B{K + ir°) x 10 6 B{K°7T+) x 10 6 



1.16±0.03±0.04 
B(K°ir°) x 10 6 



0.32 ±0.12 

F,{Ml) 



18.6 ± 1.1 

B(it + it-) x 10 6 



11.5 ± 1.3 

B(7i + n°) x 10 e 



17.9 ± 1.7 

B(tt°tt°) x 10 6 



8.9 ±2.3 



0.27 ±0.08 



5.2 ±0.6 



4.9 ± 1.1 



< 3.4 BaBar 



Table 1: Values of the input parameters used in our analysis. The CP-averaged 
branching ratios B are taken from ref. 
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Mode 


UTA 

B (1(T 6 ) 


A C p 




\V u b/V c b 
B (HT 6 ) 


A C p 




7T + 7r° 
7T 7T 


8.9 ±3.3 0.37±0.17 
5.4 ±2.1 
0.44 ±0.13 0.61 ±0.26 


8.7 ±3.6 0.39 ±0.20 
5.5 ±2.2 
0.69 ±0.27 0.45 ±0.27 


K+n° 


18.4 ±1.0 0.21 ±0.10 
10.3 ±0.9 0.22 ±0.11 
19.3 ± 1.2 0.00 ±0.00 
8.7 ±0.8 0.04 ±0.02 


18.8 ±1.0 0.21 ±0.12 
10.7 ± 1.0 0.22 ±0.13 
18.1 ± 1.5 0.00 ±0.00 
8.2 ± 1.2 0.04 ±0.03 



Table 2: Predictions for CP-averaged branching ratios B and absolute value of the 
CP asymmetries |*4cp|- The left (right) columns show results obtained using con- 
straints on the CKM parameters p and rj obtained from the UTA (the measurement 
°f The last four channels are those used for fitting the charming penguin 

parameter Pi. 



For the sake of simplicity, we also neglect here the contribution of P\ GIM . The 
B — > Kit data do not constrain this parameter very effectively, since its contribution 
is doubly Cabibbo suppressed with respect to P\. The remaining -k + -k~ mode alone 
is not sufficient to fully determine the complex parameter P^ 1 . It is interesting, 
however, to notice that the GIM-penguin contribution is potentially able to enhance 
the B{B -> 7T 7r°) up to few xlO" 6 §. 

Table ^| shows the predicted values of the CP-averaged branching ratios B and the 
absolute value of the CP-asymmetries \Acp\ for the B — > Kn and B — > nn modes, 
since the data are not able to fix the sign of asymmetries. Charming penguins are able 
to reproduce the Kn data and are also consistent with the only mr mode measured so 
far. It is interesting to notice that the latest measurements improve the consistency, 
for a comparison see refs. [§. || . 

3 Remarks on the different approaches 

Since the different approaches aiming at evaluating power-suppressed terms contain 
phenomenological parameters, it is natural to ask whether, after all, they are equiv- 
alent or not, even if the physical mechanism invoked to introduce the parameters is 
not the same. To answer this question, it is useful to compute the parameters Pi and 
pGIM w it nm improved QCD factorization. They read 

Pi = r x K a c 6 A nK + b 3 B nK , Pf /M = r x K (a c 6 - a u 6 )A nK , (3) 

where the functions aj {pi) contain the complex parameter pu (pa), see ref. for 
the definitions. These two parameters account for chirally-enhanced terms, origi- 
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nating from hard-spectator interactions and annihilations respectively, which are not 
computable within the improved QCD factorization. 

The functional dependence of the amplitudes on the phenomenological parameters 
in the two approaches is different. For instance, the GIM-penguin parameter is a pure 
short- distance correction in the improved QCD factorization, since the pn dependence 
cancels out in the difference — Og. In practice, however, the main contribution 
of the phenomenological parameters to the B — > Kir amplitudes comes from the 
annihilation term 63, i.e. from pA . This term behaves effectively as the charming- 
penguin parameter, enhancing the Cabibbo-favored amplitude. 

Notice that a vanishing pa (and pn), which turns out to be compatible with 
the data, does not mean that the phenomenological contribution is negligible. In 
fact, the parameters are defined so that the phenomenological terms are functions of 
Xa(h) — (1 + Pa(h)) ^og(m B / p h ), where the scale ph is assumed to be 0.5 GeV ||. 



4 Non-leptonic B decays and the extraction of 7 

The presence of complex phenomenological parameters in the amplitudes makes the 
extraction of 7 very problematic. We checked using the iK^/V^-constrained fit that 
almost any value of 7 is allowed, given the uncertainty on Pi. This seems a general 
problem which make us doubt recent claims proposing non-leptonic B decays as an 
effective tool for the CKM matrix determination. Even more, we think that the 
combination of the contraint from B — > Kir decays on 7 with the others can even be 
misleading. The reason is very simple: 7 is looked for through the effect of interefence 
terms in the branching ratios. The presence of a competing amplitude with a new 
phase, i.e. the one containing the phenomenological parameter, makes the extraction 
of 7 much more complicated. Although weak and strong phases can be disentangled 
in principle, in practice we checked that not only the task is very difficult now, but 
the situation improves slowly as data become more accurate, even when the CP 
asymmetries will be measured. 

Concerning various analyses based on the improved QCD factorization claiming to 
find a "large" value of 7 ~ 90°, we just notice that, as far as we know, they all assume 
the bound \pa\ < 1, suggested in ref. |J as a theoretical prejudice and supported by 
the observation that even \pa\ = produces a good fit to B(B — ► Kir). A better fit, 



however, can be obtained letting \pa\ take values up to about 3. As shown in ref. [|15 
by doing so, the contribution of the constraint from non-leptonic B decays to a global 
fit of 7 becomes totally negligible. In other words, for \pa\ ~ 3, the annihilation 
amplitude containing pa becomes competitive with the others, improving the fit to 
the Bs on the one hand and weakening the predictivity on 7 on the other. 
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